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ABSTRACT 

The origin of brown dwarfs (BDs) is still an unsolved mystery. While the standard model describes the 
formation of BDs and stars in a similar way recent data on the multiplicity properties of stars and BDs show 
them to have different binary distribution functions. Here we show that proper treatment of these uncovers a 
discontinuity of the multiplicity-corrected mass distribution in the very-low-mass star (VLMS) and BD mass 
regime. A continuous IMF can be discarded with extremely high confidence. This suggests that VLMSs 
and BDs on the one hand, and stars on the other, are two correlated but disjoint populations with different 
dynamical histories. The analysis presented here suggests that about one BD forms per five stars and that the 
BD-star binary fraction is about 2%-3% among stellar systems. 

Subject headings: binaries: general — open clusters and associations: general — stars: low-mass, brown 
dwarfs — stars: luminosity function, mass function 



1. INTRODUCTION 

Traditionally, brown dwarfs (BDs) are defined as 
(sub)stellar bodies with masses below the hydrogen burn- 
ing mass limit (HBL), niu = O.OTSM© for solar composi- 
tion, and consequent l y they cool in definitely after formation 
(iBurrows etalJll993t Ichabrier & B araffe 2000). Several at- 
tempts have been made to explain the formation of BDs by 
the same mechanisms as for stars, i. e. via fragmentation of a 
gas cloud and subsequent a ccretion (lAdams & Fatuzzoll 19961: 
iPadoan & Nordlundll2002l |2004|) . If a gas cloud in a star- 
forming region fragments there will be a certain number of 
gas clumps with masses below the HBL. Unless the mass of 
the fragment is below the local Jeans mass it will contract in 
essentially the same way as higher mass clumps and finally 
produce a single or multiple BD. This scenario predicts sim- 
ilar multiplicities and also a substellar initial mass function 
(IMF) as a continuous extension of the stellar one (e.g. the 
standard model with BDs in Kroupa et al. 2003). 

However, recent observations have shown that 
there is a lack of BD companions to lo w-mas s 

stars ([McCarthy. Zuckerman & BeckUnI l2003h . 

iGrether & Lineweaved (12006!) found a star-BD binary 
fraction among solar-type primaries of less than 1% for 
close companions, the brown dwarf desert. This implies two 
populations of stellar and substellar objects, and that binaries 
are formed in each population separately (except for pairing 
due to post-formation dynamical exchanges). Observations 
e.g. by Reid et al. (2006) also show that most BD binaries 
have a primary-to-companion mass ratio of <7 > 0.8, in 
contrast to the mass r atio distri bution of stellar bin aries which 
has typically q < 0.4 dDuquenn ov & M avodll991h. 

There are indications, e g. ' Metchev & Hillenbrandl (l2005h . 
that the BD desert may not be as dry for larger separations 
(> 30 AU) as it is for smaller ones. Since new surveys using 
adaptive optics or new instruments like the upcoming James 
Webb Space Telescope might reveal more substellar compan- 
ions to stars, the fraction of star-BD systems may increase. 

Apart from the BD desert there are more hints for 
a separate population. For example, BDs and VLMSs 
have a relatively low binary fraction of about 15% 
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parison, the stellar binary fraction is close to 100% for the 
very y oung Taurus-Auriga association (TA, about 1 Myr; 
iDuche ne 1999; Luhman et al. 2003) and about 40%-50% for 
other clusters and field stars (Lada 2006). The BD and VLMS 
binary fraction can be increased to a starlike binary fraction if 
th ere are a large fract i on of < 5 AU binaries, e.g. as deduced 
bv lJeffries & Maxtedl (1200 5). But such a semi-major axis dis- 
tribution would again imply a discontinuity of its form be- 
tween low-mass stars and VL MSs/BDs and i s not supported 
by the radial- velocity survey o fl JoergensI (12006) . We therefore 
do not consider the starlike formation as a major mechanism 
for BDs. 

It has also been argued that the low binary fraction of 
BDs can be understood as a continuous extension of a trend 
that can already be recognized from G d warfs to M dwarfs 
(lLuhmanll2004al: ISterzik & Durisenll2003h . Therefore, the bi- 
nary fraction alone cannot be taken as a strong evidence to 
introduce a separate population. 

The most striking evidence for two separate populations 
is the empirical fact that the distribution of the separations 
and therefore the binding energies of BD bin aries differs sig- 
nificantly from that in the stellar regime (Bouv et aL '2003j; 
Burgasser et al. 2003; Martin etal. 2003; Close et al. 2003). 
This is shown in Figures [T] and |2l The solid line is a 
Gaussian fit to the central peak of the histogram while the 
dash-do tted one refers to dBasri & ReinersI 120061 a com- 
pressed |Fischer&^arc3[i993^t)7^Dr'an^ very low-mass 
stars (VLMS) have a semi-major axis distribution limited 
to < 15 AU, whereas M, K, and G dwarfs have a very 
broad a nd similar distributi o n (long-dashed and short-dashed 
curves: IClose et al.l|2003l: iLaw. Hodgkin & Mackavl 120071; 
iGoodwin et al.l I2007IK There is also a dearth of BDs be- 
low 1 AU. Recent findings, e.g . by iGuenther & Wuchterll 
( l2003h and iKenvon et alj d2005h . suggest a low number of 
such very close BD/VLMSs binaries. The semi-major axis 
distribution ofBDsA'^LMSs binaries based on the data from 
IClose et all d2003h can be modelled with a log a Gaussian 
centered at 4.6 AU (loga = 0.66) with a half-peak width 
of cr = 0.4. It corresponds to an overall BDA'LMS binary 
fraction o f /bd- bd = 0-15. If data fro r n , Lu hman (2004bi) , 
iJoergensI d2006l) . and iKonopackv et alj d2007h are taken as 
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hints to inc omplete data between ab out 0.02 and 1 AU, the 
co mpressedlFischer & Marcyl (Il992h Gaussian from Figure 4 
in iBasri & ReinersI (l2006l) may provide an appropriate enve- 
lope. However, for an assumed BD mass of 0.07 Mq and 
/bd-bd = 0.26 their period distribution corresponds to a semi- 
major axis distribution with a w 0.85 and is therefore still in- 
consistent with that of M and G dwa rfs. 

Although lKonopackv et al.l (|2007|) have recently found five 
VLMS binaries in TA with four of them having separations 
much larger than 15 AU, the sudden change of the orbital 
properties remains. In particular, they found two binaries with 
projected separations slightly above 30 AU and two others 
with separations between 80 and 90 AU. P ossib le implica- 
tions of these discoveries are discussed in § 15.41 However, 
the truncation near 15 AU cannot be derived from the stellar 
distribution through downsizi ng according to Newton's laws 
dClose et al . 2003; B ouv et al.ll2003l) . Not even dynamical en- 
counters in dense stellar environments can invoke such a trun- 
cation n ear 15 A U (Burgasser et al. 2003). Using A^-body sim- 
ulations Kroup a et al.l (120031) tested the hypothesis that BDs 
and stars form alike, and mix in pairs, and found that, de- 
spite of close dynamical encounters, the distribution of the 
semimajor axes of BD binaries remains starlike. That is, dy- 
namical encounters even in dense clusters cannot truncate the 
BD binary distribution near 15 AU. They further found that 
star-BD binaries would be much more frequent than actually 
observed. Thus, this classical hypothesis is rejected with high 
confidence. 

With this contribution we study the implications of the ob- 
served change of binary properties on the underlying single- 
object initial mass function (IMF). For this purpose, we an- 
alyzed the observational mass functio ns of TA (Luhman e t 
al. 2003a, 2004b), the IC 348 cluster (iLuhman et alJ lIOOl). 
the Trapezium cluster dMuench et al.ll2002l). a n d the Pleiades 
cluster based on data by i Dobbie et alj (|200"2|) . |Moraux et al.' 
(|2003), and the Prosser and Stauffer Open Cluster Database.^ 
For all the systems we analyst, we refer to the MF as the IMF, 
although strictly speaking this is not correct for the 130 Myr 
old Pleiades. In §|2]we shortly review the definition of the 
IMF and the role of multiplicity on the shape of the individual 
body IMF compared to the observed system IMF (IMFsy,). 
We show that the IMF that can be derived from an observed 
IMFsys does not need to be continuous even if the IMFsy, does 
not show any discontinuity. In § |3] we describe the statisti- 
cal method how to fit a model IMFsys by combining BD and 
star IMFs as an approximation to the observed IMFsys. § E] 
presents the results that indicate a discontinuity close to the 
HBL. Also the BD to star ratio from the model is calculated 
there. In §|5]the results are discussed in the context of four 
alternative BD formation scenarios, i.e. embryo ejection, disk 
fragmentation, photoevaporation, and ejection by close stellar 
encounters. 

2. THE IMF FOR INDIVIDUAL STARS AND SYSTEMS 

2.1. Definition 

The IMF is among the most important properties of a stellar 
population since it gives hints to the processes that form stars 
and BDs. Although we can only observe stellar populations 
at their given age there are data of several very young popu- 
lations where the mass function is probably still very close to 
the initial one. 

Available at |http://www.cfa.harvard.edu/~stauffer/openci7] 
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Fig. 1. — The semi-major axis {a) distribution for VLM bi- 
naries (VLMBs) with a total mass m[,in < O.IMq based on data 
from Nick Siegler's Vei-y Low Mass Binaries Archive (VLMBA, 
http://paperclip.as.arizona.edu/~nsieglerA'LM_binaries/, version from June 
1, 2007; histogram). The data can be fitted with a log-normal distribution. 
For the VLMBs (solid curve) the Gaussian parameters are ,tpcak = 0.66 (cor- 
responding to fl = 4.6AU) and cr = 0.4 witli a normalization factor (= total 
binary fraction) /tot = 0.15, the binary fraction among BDs and VLMSs. For 
the M0-M4 dwarf binaries by Close et al. (2003) ipeak = 1-6 and = 1.26 
(/tot = 0.3) and for the G dwarf binaries Xpe^j^ = 1.53 and cr = 1.53 (/tot = 
0.5; Fischer & Marcv 1992, long-dashed and short-dashed curves, respec- 
tively). For VLMBs the gap within -1.5 < \oga/AU < can be interpo- 
lated with the compressed Gaussian of Basri & Reiners 1 2006) with cr = 0.85 
(thin dot -dashed curve), also enveloping wider BDA'LMS binaries found by 
IKonopac kv et al. (2007). It corresponds to a total BDA'LMS /bd-bd = 0.26. 
The areas below the histogram and the curves are equal to the conesponding 
/tot (eq.|3). 
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Fig. 2. — The distribution of the orbital binding energy iibind = 
-0.5G(m\m2l a) for VLM binaries and low-mass stellar binaries. Line types 
and symbols are the same as in Fig. [T] As for the semi-major axes there is a 
clearly different distribution for BDsA'LMSs on the one hand side and M, K 
and G dwarf binaries on the other 

In general, the mass distribution of stars and BDs can be ap- 
proxi mated by a powe r law or a combination of several power 
laws. ISalpeted d 19551) estimated the relationship between the 
stellar mass and the relative number of stars of a given mass 
as a single power law. 



An 

t,{m) = — = km 
dm 



(1) 



for 0.4 < hi/Mq < 10 and with a = 2.35 and a normalization 
constant k. The IMF is often expressed on the logarithmic 
mass scale and then becomes 



^L(logm) : 



dn 



dlogm 



: (In 10) m^(m) = ki^m 



(2) 
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where k]^ is the corresponding normahzation constant for ^l- 
An IMF is said to he flat if ^L(logm) is constant, i.e. a= I. 
More recent work has shown there to be a flattening in the 
lower-mass regime of the observed mass function and, in the 
representatio n, even a turnover nea r the BD-star transition 
kroupa 2001: Reid. Gizis & Hawleyll200 2: Ch abrier 200 3). 
Here, all objects from BDs to the most massive stars in a 
cluster are described by a continuous IMF, i.e. with a single 
population containing BDs as well as stars. The universal or 
"canonical" IMF has ago = = 0.3 (0.01 < m/Mg < 0.075), 
ai = 1.3 (0.075 < ot/Mq < 0.5), and ^2 =2.3 (0.5 < m/MQ < 
nimux), where m,nax is giv en by the mass of the host cluster 
(IWeidner & KroupalllOO^ . The generally accepted wisdom 
has been t hat the IMF is continuous from above 0.01 Mq to 
m^ax re.g. lChabriei<l2002h . 

2.2. Unresolved Binaries 

Star cluster surveys are usually performed with wide-field 
telescopes with limited resolution that do not resolve most of 
the binaries. Hence, they yield, as an approximation, system 
MFs in which unresolved binary systems are counted as one 
object. The fraction of unresolved multiples can be taken as 
the total binary fraction of a cluster. 



/to 



(3) 



where A^sng is the number of singles (or resolved individual 
bodies) and A'bin.tot the number of (unresolved) binaries. Un- 
resolved binaries increase the number of individual bodies, 
N\,od, in a cluster such that 



A^bod = (l+/tot)M 



sys ; 



(4) 

where A^sys = Nsng +A^bin.tot is the total number of systems. Note 
that a "system" is either a single body or a binary or a higher 
order multiple. Her e we ignore higher order m ultiples, be- 
cause they are rare dGoodwin & Kroupall2005l) . The binary 
fraction in dependence of the mass of the primary star, mpnm, 
is 

\ ^bin,tot(^prim) 

/(Wprim) = T^— —r: z r , (5) 

A'sng(Wpi-ini) + A'bin,tot(Wprim) 

where A^bin.tot('«prim) and A^sng('«prini) are, respectively, the 
number of binaries with primary star mass mprim and single 
stars of mass mpi-im. 

2.3. The System Mass Function 

The effect of this binary error on the appearance of the IMF 
can be described as follows. Assume a stellar population with 
/tot < 1 and with stellar masses with a minimum mass m^i„ 
and a maximum mass limit of m^^^. The minimum mass of 
a binary is 2mj„i„, while the mass-dependent binary fraction 
/(m) = for mmin < 'Wsys < 2OTjnin- A binary closely above 
2minin can only consist of two stars near m„^a making such 
binaries rare. For higher system masses, where a system can 
be a multiple or a single star, there are more possible com- 
binations of primary and companion mass, so that the binary 
fraction increases with the system mass and approaches an 
upper limit for the most massive objects. 

Figure |3] shows the general shape of a system IMF for a flat 
(logarithmic scale) IMF with m^i„ = Q.IMq, m„iax = ^Mq, 
and /tot = 0.5. The IMFjys is flat between m^i„ and 2mjai„ 
(logm/M© « -0.7) and rises above 2m,^in to a maximum at 
'«max- Systems with Wsys < 2m„iin can only be singles and the 
IMFsys in this region is just the IMF minus the mass function 



of objects that are bound to a multiple system. For masses 
msys > OT,nax, on the other hand, only binaries exist, and the 
IMFsys declines towards zero at mjys = 2m,nax, the highest 
mass possible for binaries. The sharp truncation of the IMF at 
m = mmax causes the sudden drop at msys = Wmax, while the mi- 
nor peak at msys = '«max + Wmin corresponds to the maximum 
of the binary mass function (IMFbm, the IMF of binary sys- 
tem masses). For a binary of this mass the primary and com- 
panion mass can be drawn from the whole IMF, and thus the 
number of possible combinations becomes maximal. It should 
be noted that natural distributions with smoother boundaries 
probably do not show such a double peak. 

Mathematically and in the case of random 
pairing, which is a r easonable approxim atio n in 
the stellar regime (jM alkov & Zinneckej 120011: 
iGoodwin. Kroupa. Goodman & B urkert 2Q(yT), the bi- 
nary mass function is just the integral of the product of the 
normalized IMF, ^ = ^/A^bod, of each component times the 
total number of binaries A^bin.tot- Given the masses of the 
binary components A and B, nip, and m^, the binary mass 
fnbin = m\ + mB. Thus, m^ = mbm — m^. Thus, IMFbin can now 
be written as 



Cbin(»Ibin) = M 



bin.tot 



^(m)^(mbi„-m)dm, (6) 



where mmn is the lower mass limit of all individual bodies in 
the population. The upper limit of the integral, mbi,, - m^in, 
is the maximum mass of the primary component m\ corre- 
sponding to a secondary component with m^ = Wmin- 

The other extreme case of assigning the component masses 
is equal-mass pairing. In that case, equation ^ simplifies to 
the IMF of one of the components and IMFbi,, is just the IMF 
shifted by a factor of 2 in mass and corrected for binarity using 
equations ^ and (|4|i: 



/tc 



^binequal('«bin = 2m) = ^(m) 
t + ./tot 



(7) 



This case is more applicable for BD binari es since the ir mass 
ratio distribution peaks at a ratio q = I dReid et alj |2006). 
However, due to the low overall binary fraction of BDs the 
effect of the mass ratio distribution is quite small. 

The IMFsys is just the sum of IMFbin and the IMF of the 
remaining single objects, (1 -ftodS,- 

£,sys(m) = 6in(ni) + (l - /tot) ^(m) . (8) 

Thus, to obtain the true individual star IMF, ^(m) has to 
be extracted from the observed ^sys('w) for which a model for 
Cbin('«) is required. That is IMFsys has to be corrected for 
unresolved binaries in the cluster or population under study. 
This leads to a significant increase of the numbers of ob- 
jects at the low-mass end because low-mass objects contribute 
to both low-mass singles and intermediate-mass binaries and 
thus more individual objects are required to reproduce the 
observed IMFsys (Kroupa etal. 1991; Malkov & Zinneckej 
1200 lb . In the mass range mi„in-2mmi„, systems can only be 
single because the system mass is the sum of the masses of 
the system members. Increasing the system mass beyond 
2mmin causes the binary contribution to rise quickly and then 
to asymptotically approximate a maximum value. Thus, the 
fraction of singles among M dwarfs is higher than for G 
dwarfs (tLadaL2006,) . as one would expect for a stellar pop- 
ulation (> O.O8M0), i.e. essentially without BDs. 
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Fig . 3 . — The system IMF (IMFsy.s solid line) corresponding to a flat loga- 
rithmic IMF (i.e. a = 1 , dashed line) for M7min = 0.1 <m/MQ < 1 — "^nia.\ with 
an overall binary fraction of 50%. The con'esponding binary IMF (IMF(,i„ eq. 
|6) is shown by the thin dotted line. The system MF peaks at 1 Mq and is trun- 
cated for higher masses while there is a minor peak at 1 . 1 Mq corresponding 
to the peak of the binary MF at m^^ + mmax. 
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Fig. 4. — A population of stars and BDs with different binary fractions 
can result in a discontinuous IMF, even if the observed (system) IMF appears 
to be continuous. The binary fraction of the BDs is lower than that of stars 
and therefore a lower number, Wbod. of individual objects is required for the 
frequency of BD and stellar systems of a given mass being equal (eq. |4). 
Note the overlap region indicated by the horizontal dotted lines: some starlike 
bodies may actually be physical BDs (upper line) while some BD-likes are 
indeed very-low-mass stars. 

Figure |4] schematically shows the effects of correcting 
a flat continuous observed IMFjys for unresolved multiples 
whereby the binary fraction for BDs is smaller than for stars, 
as is observed to be the case. For a binary fraction of 
15% (i.e. /bd-bd = 0.15) among BDs there would thus be 
A^BD = (l+/BD-BD)A^sys,BD = l-15A^sys,BD BDs in total, whilc 
for stars with /tot = 0.5 we would have 1 .5A^sys.stai- This exem- 
plifies how the change in / leads to a discontinuous IMFsyj. 

3. THE METHOD 

3.1. The Parameter Space 

The most straightforward way to calculate the influence of 
binaries on the s tellar/substellar statistics would be the Monte 
Carlo method ( Krou pa et all 1 1 99 iL [1993 1) . However, for bet- 
ter (smoother) results and to reduce the computational efforts 
we do not use a Monte Carlo approach here but a semiana- 
lytical approach in which the binary mass function is calcu- 
lated via numerical integration of eq. ^ for each population. 



TABLE 1 

List of the variable parameters that span the parameter 

SPACE OF a two-population IMF STUDIED IN THIS CONTRIBUTION. 

Note that qbd is not varied for IC 348 and the Pleiades since 

IT IS NOT well-constrained BY THE DATA FOR THESE CLUSTERS. 



Parameter No. 


Quantity 


varied from 


to (approx.) 


1 


"BD 





1 


2 


'"max.BD 


n«o,stiu- 


'no.stiu+0.2 


3 


log7?.pop 


-1 






TABLE 2 

Size of the observed sample, Wsys (= the number of observed 

SYSTEMS), NUMBER OF DATA POINTS (MASS BINS), njaia, NUMBER OF 
FITTING PARAMETERS, C, AND THE NUMBER OF DEGREES OF FREEDOM, 
I' = "data - 1 -C, FOR THE FOUR CLUSTERS UNDER STUDY. 



Cluster 




"data 


c 


V 


Trapezium 


1040 


30 


3 


26 


TA 


127 


9 


3 


5 


IC 348 


194 


9 


2 


6 


Pleiades 


~500 


19 


2 


16 



This can be done with a standard quadrature algorithm. This 
algorithm has been verified with a Monte Carlo simulation 
of a few million random experiments, each being a random 
draw from the IMF. The Monte Carlo method is used later to 
determine the BD-to-star ratio and the total binary fractions 
within defined mass ranges (§§ 14.21 and 14.31 ). It makes use 
of the Mersenne Twister random num ber generator developed 
by iMat sumoto & Nishimural (Il998h . Since there are only a 
few runs to be done (in contrast to hundreds of thousands of 
runs during an iterated parameter scan, see below) the sim- 
ple Monte Carlo approach with an appropriately large random 
sample is fully sufficient for this purpose. 

The lack of star-BD binaries and the truncation of the BD 
binary separation distribution suggest two disjunct popula- 
tions where binary components are taken only from the same 
population. For reasons discussed in §|5]we call these popu- 
lations BD-like and starlike. The binary corrections are there- 
fore applied separately to the stellar and the substellar regime 
rather than to a combined population. Random pairing over 
the whole mass regime is n ot considered furthe r as it leads to 
too many star-BD binaries dKroupa et al.ll2003l) . 

The approach here requires separate application of normal- 
izations on both populations. For this purpose we define the 
population ratio, 

np,p=Ni/N2 (9) 

where A^i is the number of individual bodies of the BD-like 
population (wjo.bd <m< Wmax.Bo) and A^2 ('"o.star <m< m^ax) 
that of the starlike population. This must not be mixed up with 
the BD-to-star ratio, TZ, (eq. lfT4l ) which refers to physical 
BDs and stars separated by the HBL. 

The partial IMF for each population can be described by 
parameters a,, nij, and a normalization constant. In this work 
a single power law for BDs (/ = "BD") and a two-part power 
law (;=1 or 2) for stars is applied. Thus, there is a mass border, 
mi2, separating the two power law regimes. The parameters 
of both populations form a three-dimensional parameter space 
of the IMF model for each cluster. It has been found that the 
lower mass limits of the BD-like population, mo BD = 0.01 M©, 
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and that of the starHke population, mo.star = O.OTM©, are 
suitable for all studied clu sters. Further more, we focus on 
the canonical stellar IMF (lKroupall2()0lb with ai = 1.3 for 
m < Q.SMq = mi2 and = 2.3 for m > O.SM©. The reason is 
that this canonical IMF has been verified with high confidence 
by other observations as well as theoretically. Thus, only the 
BD-like power law, aBo, the normalization of the BD-like 
population against the stellar one, T^pop, and the upper mass 
border, Wmax.BD^ of the BD-like regime are the parameters to 
be varied. Because of the sparse and probably incomplete 
data sets for IC 348 and the Pleiades, qbd has also been set 
to the canonical value (i.e. aeo = ao = 0.3) for these clusters, 
varying only 7^pop and mmax.BD- 

Table [T] lists the variables and their range of variation for 
the Trapezium, IC 348, and the Pleiades. As the upper mass 
limits of stars in the clusters one can either take the max- 
im um observed mass o r a th eoretical mass limit, as given 
by IWeidner & Kroupal (l2006b . Because in our model the 
IMF is cut sharply rather than declining softly as shown in 
that work, we here set the mass limit somewhat below the 
IWeidner &^ Kroupa (2006) limits. For the Trapezium we set 
an upper limit of 10 Mq, 1.5 Mq for TA, 3 Mq for IC 348 
and 5 M© for the Pleiades. Note that the observed most mas- 
sive star in the Trapezium has a mass of w 50M©, but varying 
'Wmax to this valuc does not affect the results significantly. 

Little is known about nonplanetary substellar objects below 
the standard opacity limit for fragmentation, 0.0 IM© , mak- 
ing this a reasonable lower mass limit for our BD statistics. 
The lower-mass limit for stars, however, is chosen somewhat 
arbitrarily but is justified by a smaller test study that finds less 
agreement with observational data used here if the stellar mass 
is lowered too much below 0.08Mq. As will be shown later, 
the extreme case of a stellar mass range that includes the BD 
mass range leads to poor fits to the observational MFs as well 
as to some inconsistency with the observed binary fraction. 
This is discussed later in §§|4]and|5] In general, the choice 
of the mass ranges also influences the mass ratio distribution, 
since for random-pairing each component is distributed via 
the IMF. Thus, for a given primary mass out of a chosen pop- 
ulation, the companion mass distribution is just equal to ^(m) 
within the mass interval m^in < m < mp^n,. 

The overall binary fraction is taken as a constant for each 
population. For the stars we adopt an unresolved binary frac- 
tion of 80% for TA (/star-star = 0.80) and 40% for the others 
(/stai-stai = 0.40), in accordance to the observational data. For 
BDs we choose a value of 15% for all cases (/bd-bd = 0.15). 
Furthermore, we found the need to introduce an overlap of 
the BD and star regimes between 0.07 and about 0. 15-0.2Mq. 
Indeed, there is no reason why the upper mass border of the 
BD-like population should coincide with the lower mass bor- 
der of the starlike population. The physical implications of 
this required overlap region are discussed in § |5] 

3.2. Minimization 

In § 13.11 a parameter space of three (Trapezium, TA) and 
two (IC 348, Pleiades) dimensions of individually adjustable 
quantities has been defined. For each cluster, the quality of a 
set of fitting parameters is characterized via the x' criterion, 

x' = v'"'-"V""'>\ (.0, 

I ' 

where 'Csys('Wi) is IMFjys at the midpoint of the ith bin of the 
observational histogram that is taken as the reference and «data 



is the number of mass bins in the observational histogram 
(Tab. |2]i. The values cr, represent the error bars of the ob- 
servational data, or, if none are given, the Poisson errors, and 
A', is the number of systems found in the ith mass bin. The 
fitted IMFsys is normalized to the total number of systems in 
the cluster. 

That set of c fitting parameters that minimizes the reduced 
X^ value, xt = X^/^^^ defines our best-fit model. Here = 
«data~c - 1 is the number of degrees of freedom. The values 
for Hdata, c, ly, and the sample sizes for the studied clusters 
are listed in Table |2] The reduction of by 1 is due to the 
normalization of the total number of stars and BDs against the 
observational data. For the two-component model the fitting 
parameters are the power in the BD regime, aso, the upper 
mass border of the BD population, and the ratio TZpop', thus, 
c = 3 (2, for IC 348 and the Pleiades). The probabiUty, P, 
that the model has a reduced x^, xi = X^/'^' large as or 
larger than the value actually obtained is calcu lated from the 
incomplete Gamma function (iPress et alJll992h . 

Q(xlW) = P(xl<xl)- (11) 

The logarithmic mass error is not mentioned in the sources, 
but is given by the photometric measurements and to a larger 
degree by uncertain theoretical models of stars and BDs. We 
assumed a value of Alogm = 0.05 for the Trapezium and 
Alogm = 0.1 for the others, corresponding to a relative error 
in the mass estimates of 12% and 26%, respectively. This er- 
ror has been estimated from the width of the logarithmic bins 
in the observational data for the Trapezium, TA and IC 348. 
For the Pleiades, for which non-equally spaced data from dif- 
ferent sources are given, we assume the same log mass error 
as for TA, and IC 348. To take this into account the logm 
values have been smoothed by a Gaussian convolution corre- 
sponding to a lognormal smearing of the masses. 

After generating the model IMFs for both populations a bi- 
nary mass function (eq. ID) is derived separately from the 
stellar and substellar IMF, i.e. there are formally no BD-like 
companions to starlike primaries, in accordance with the ob- 
servations. This leads to consistency with the observed binary 
fraction (§ 14. 3b and the BD desert. Note, though, that as a re- 
sult of the required overlap both the BD-like and the starlike 
population contain stars as well as BDs; thus, we do have star- 
BD pairs in our description (more on this in §|5]). 

The addition of the two resulting system IMFs leads to 
an overall IMFsys for the whole mass range. By adjusting 
the power law coefficients and the population ratio 7?^pop, the 
IMFsys is fitted against the observational data such that x^ is 
minimized. The promin ent substellar peak i n the Trapezium 
cluster below 0.03 M© dMuench et alj|2002l) is ignored here 
because it is well below the BD-star mass limit and there- 
fore does not interfere with any feature there. Furthermore, 
it is possibly an arte fact of the BD mass-luminosity relation 
(iLada & Ladall2003h . 

3.3. Error Estimation 

The errors of the parameters are estimated from the 
marginal distribution of each parameter within the parame- 
ter space (lLuptonlll993h . The marginal probability density 
distribution, p, is 

p(Ai) = dP(Ad/dA, 

= 11 £(AuA2,Ai)dAkdAj, ^'^> 
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Fig. 5. — The marginal probability density distributions of the parameters 
obd, logT^pop and log (ni max.B p/^n ) as an estimate of the errors for the 
Trapezium cluster fit (see § 13.31 . The peak-width engulfing 68% of the whole 
parameter sample (thus referring to Ki Icr, dashed vertical lines) is taken as 
the error for each parameter Because the median (solid vertical line) of the 
sample does not match the best-fit value (dash-dotted line) in all cases, both 
are listed in Table[3] 



where Ai,A2 and A3 correspond to q;bd, Tlpop, Wmax.BD and 
i= 1, ... ,3; j,k^i. The HkeHhood C is defined as 

C{AuA2,Ai) = Me-^-l\ (13) 

where the normahzation constant M is such that the integral 
over the whole range is one. Note that instead of is 
used here. The integral is approximated by summation with 
equidistant stepping of each parameter We conservatively re- 
jected all parameter sets with P < 0.27% (corresponding to 

As an example, the parameter distribution for the Trapez- 
ium is illustrated in Figure |5] The interval around the median 
which contains 68% of the scanned parameter sets is taken as 
the a measure of the errors. The interval limits and the me- 
dians are shown among the best-fit results in Table |3] Appar- 
ently, the median of the probability density distribution does 
not always coincide with the best-fit value. The reason for 
this is that the distribution is slightly asymmetric for most pa- 
rameters and clusters with the Pleiades being the by far worst 
case. This results in asymmetric error bars. Note that there 
are sets of a^D, ^pop, and mmax.BD that are within these error 
limits but with a > 3cr, and thus the error limits may be 
slightly over-estimated. 

Another possibility to illustrate the statistical significance 
of a fit is by its confidence contours within a two-dimensional 
subspace of the parameter space, as is shown later in § 14.21 




-2 -1.5 -1 -0.5 0.5 
log-|o m//W^, 



Fig. 6. — Four different IMF models applied to the Trapezium cluster. Top 
panel: A classical continuous-IMF fit (dashed line) for the Trapezium cluster 
with random pairing over the whole BD and stellar mass range. The method 
of converting the IMF to the IMFsys (solid curve) is essentially the same as 
for the two-component IMFs discussed in this paper but with no allowance 
for an overlap or a discontinuity in the BD-star transition region. Although 
the general shape of the observational histogram (thin steps) is represented 
by this fit this model leads to an unrealistic high fraction of star-BD binaries 
and is therefore discarded. Second panel: The Trapezium IMF fitted with 
two separate IMFs, §bd (dotted fine) and ^siar (dashed line), as in Figure |7] 
(the two-component IMF). Third panel: The same but with the sum IMF of 
both component IMFs shown here as the possible appearance of the IMF if 
all binaries could be resolved. Apparently, there is a "hump" between about 
0.07 Mq and 0.2 Mq bracketing the probable overlap region of the BD-like 
and starlike populations. Bottom panel: Trapezium fitted with two separate 
IMFs but with no overlap. The binary correction leads to the dip near the 
hydrogen-burning limit. However, the fit is only slightly worse than that one 
with an overlap. 

4. RESULTS 

4. 1 . The IMF for BDs and Stars 

Figure|6]demonstrates the results for different models (con- 
tinuous IMF, two-component IMF with and without over- 
lap region) for the Trapezium cluster. The continuous IMF 
is shown for illustration only. Its underlying model as- 
sumes a single population containing BDs as well as stars 
(Kroupa et al. 2003). Instead of an overlap as for the two- 
component model the mass border between the regime of the 
BD slope Q!BD and the (canonical) stellar regime is varied. The 
total binary fraction, /tot, is set to 0.4 with random pairing 
among the entire population. Our calculations give a best fit 
with aBD = -0.4 ± 0.2, m„,ax.BD = '«o,star = 0.093 ± 0.01 (and 
canonical stellar IMF) which is comparable to the canoni- 
cal IMF for BDs and stars (Kroupa 2001) . As already men- 
tioned in § [H it leads to a high number of star-BD binaries. 
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TABLE 3 

The best-fit BD-like power law coefficients, a^u, population ratios, 7?.pop, and BD-like upper mass limits, m„,ax,BD- The 

UNCERTAINTIES ARE DERIVED IN §[3] NOTE THAT ObD IS SET TO THE CANONICAL VALUE FOR IC 348 AND THE PLEIADES. 



Cluster 


Median ±la 

QBD log Tlpop nimax.BD 


Best fit 

QBD log7?-pop "Jraax.BD 


xi P 


Trapezium 
TA 

IC 348 
Pleiades 


0.47^:1^ -0.67:0;i0 0.22^:0* 
-().08!?;g -0.68!«;|3 O.U^f, 
- -0.57:«;?3 0.22^06 
-0.82!«:tt 0.07!»i 


0.45 -0.66 0.22 
0.00 -0.62 0.12 
0.3(c) -0.51 0.22 
0.3(c) -0.80 0.06 


0.23 0.99998 
2.69 0.020 
1.01 0.416 
1.20 0.257 




ape 

ters based on observational data by 'Luhman! (2004a!), 'Luhman et al.' (20031) 
and Muench et al. (2002) (solid histograms), as well as for the Pleiades clus- 
ter using data from M oraux et al . 12003) (filled circles), Dobbie et al . ( 2002) 
(open squares) and from the Prosser and Stauffer Open Cluster Database 
(from which only those for logm/MQ < 0.55 are used for fitting). For all 
clusters the stellar IMF, shown here a s the dashed line, has the canonical 
or standard shape after KrouplJ 120011) , i.e. Oj = 1.3 below 0.5 Mq and 
02 = 2.3 above 0.5 Mq. The Trapezium IMF is fitted without the substellar 
peak indicated by the faint dotted steps. The dotted lines show the optimized 
binary-corrected BD-Uke body IMFs (distribution of individual BD-like bod- 
ies). The resulting model system IMFs are represe nted by t he solid curves. 
Ages are from Luhman et al. 1 2003), Luhma n efaTl )2003() , lHillenbrand et alj 
(^01 ) and Barrado y Navascues et al. 1 2004) for TA, IC 348, Trapezium and 
the Pleiades, respectively. Note that the BD-like IMF of IC 348 and the 
Pleiades are set to the canonical slope qbd = 0.3 since they are not well 
constrained due to the lack of BD data for these clusters. 



O 
T3 

•D 



O 




P = 0.23 



-1 -0.5 
logio m/Mr, 



0.5 



Fig. 8. — The two-component model with canonical stellar IMF for a 
slightly modified TA observational histogram (black soHd steps). The peak 
at 1 Mq and the dent around 0.3 Mq have been slightly flattened without 
leaving the error bars of the original data (grey steps and bars). This fit has 
a reasonable confidence of about P = 0.25, about 12 times larger than for the 
original histogram (P = 0.02, see Tab. [3j. 



/stai-BD = 15% ±4% altogether (see j43b . while lKroupa et alj 
0)03) found 7%-15% star-BD binaries with a > 30 AU in 
their standard model with stars and BDs. Furthermore, the 
different orbital properties of BD-BD and star-star binaries 
(Fig. [1) cannot be explained by a single-population model 
without further assumptions. Therefore, it is not considered 
further here. 

The two-component IMF model (Fig. |6] second panel) ac- 
counts for the empirical binary properties of BDs and stars. 
Apart from this it also fits the observational data (Fig. |6] his- 
togram) slightly better than the continuous model, especially 
in the BDA'^LMS region and the "plateau" between 0. 1 and 

0.5 Mq. 

Objects with equal mass and composition appear as equal 
in observations, even if they have been formed in different 
populations. Thus, it is hard to determine their formation his- 
tory. A high-resolution survey that resolves most of the mul- 
tiples would yield an overall mass function from the lowest 
mass BDs to the highest mass stars. However, if such a re- 
solved individual mass function is composed of two overlap- 
ping populations, it may be possible to detect its imprints as 
an excess of objects within the overlap region. To illustrate 
this ((Fig. |6] third panel), we construct such an overall IMF 
by simply adding the two IMF components from the second 
panel. BDs and stars are paired separately from their respec- 
tive mass range. It is clear that the addition of the two IMFs 
leads to discontinuities, i.e. in the overlap region the contin- 
uous IMF increases steeply at the minimum stellar mass and 
drops again at the upper mass limit of the BD-like regime. 
Even with a smoother drop at each end of the BD-like and 
starlike IMF it is likely to be detectable once most of the 
multiples have been resolved. Thus, a discovery of such a 
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Fig. 9. — Contour plots of 5% and 1% significance levels from fitting 
(see §|3]for details) of obd ™d the BD-like to starlike ratio 7?.hbl at m = 
0.075 Mq (eq. [B) for the Trapezium, TA, IC 348 and the Pleiades. The 
stellar IMF is canonical, i.e. ai = 1.3 and 02 = 2.3 while the upper BD mass 
limits, mma,K,BD. are the best-fit values from Table [3] All fitting parameters 
outside the solid line are rejected with 95%, and those outside the dashed line 
are rejected with 99% confidence. The optimum of the Trapezium and TA 
is marked by the filled circle while the cross mai'ks the standard/canonical 
configuration with a continuous IMF (i.e. log7?.HBL = 0) and cjbd = 0.3. As 
can be seen, it is well outside both levels (even for arbitrary cvbd)- For the 
Pleiades it is still well outside at least for reasonable choices of c^bd ■ Note 
that the optimum for obd for IC 348 and the Pleiades does exist in this 2D 
cross section but not in the full 3D parameter space. 



0.8 
0.6 
0.4 
0.2 

1 

0.8 
0.6 
0.4 
0.2 

1 

0.8 
0.6 
0.4 
0.2 

1 

0.8 



Taurus 
1 IVIyr 




IC 348 
2 Myr 




Trapezium 
1.5 Myr 




Pleiades 
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Fig. 10. — The binary fraction / (eq.|5) of the four discussed clusters as a 
function of primary mass, niprim, for the two-component best-fit IMF model 
(with canonical stellar IMF) is shown by the solid line. The BD binary frac- 
tion appears flat due to the equal-mass pairing used in our algorithm. Since 
merely very sparse data is available for low-mass BDs, only the higher-mass 
end of the curve should be taken into account. For comparison, the average 
binary fractions of G, K, M and L type stars in the solar neighbourhood (av- 
erage age ~ 5 Gyr; filled circles: Lada 2006, open squares: Kroupa et al] 
2003) have been added to each panel. The shaded region in the top panel 
indicates the approx. 1 Myr old pre-main sequence data (Duchene 199^ 
while the thin solid line represents the initial TA-like dynamical model from 
Kroupa et al. (2003). The single fill ed square in the top panel at log m^-1.2 
is the VLMS/BD datum inferred bv lKraus et al. ( 2006) for TA. Note that the 
stars in TA have a binary fraction near 100% tDuchene 1999. : Kroupa et afl 



"hump" within the resolved IMF would strongly support the 
two-component model. 

To illustrate the effect of the overlap region, a two- 
component IMF without an overlap is shown in the bottom 
panel of Figure |6] Like the continuous IMF this approach 
gives a slightly worse fit, especially between 0. 1 and 0.2 Mq 
(-1 < logm/M0 < -0.7). However, the "dip" in this region 
is only weak. 

The results of the two-component IMF model for all clus- 
ters are shown in Figure |7] the power law coefficients being 
hsted in Table [3] For the Pleiades and IC 348 the BD IMF 
slope has been kept standard, i.e. qbd = 0.3 (Kroupa 2001), 
because the sparse, and in the case of the Pleiades most prob- 
ably incomplete data, do not allow useful confidence limits to 
be placed on qbd- For the same reason the highly uncertain 



observational data for log(m/M0) > 0.65 have not been used 
for fitting (but are still plotted for completeness). 

As can be seen the BD/VLMS and stellar IMFs do not meet 
at the BD-star boundary. The number density of individual 
BDs near the BD-star border is about one third of that of indi- 
vidual low-mass M dwarfs. This discontinuity cannot be seen 
directly in the observational data because it is masked by the 
different binary fractions for different masses ( ^2.3l l. but the 
IC 348 data (histogram in the third panel of Figure |7]i may 
indeed be showing the discontinuity (compare third panel of 
Figure |6]l. 

Furthermore, there are large uncertainties in the mass de- 
termination of stars in stellar groups such as TA which may 
probably lead to even larger error bars of the bins than shown 
in the observational data. Even within the given error bars 
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certain variations of these data provide a much better fit to 
the canonical stellai- IMF with P = 0.25 (original: P = 0.02, 
Table O. This has been done in Figure |8] by lowering the 
peak near 1 Mq and somewhat rising the "valley" around 
0.3 M0 corresponding to the re-shuffling of 9 stars out of 127 
(Tab. [TJ by 0.1 to 0.4 M© (one bin width), e.g. through 
measurement errors . This suggests that, in agreement with 
iKroupa et al.l (l2003h . the TA IMF might not necessarily be in- 
consistent with the canonical IMF. 

Except for the Pleiades all clusters show features near the 
peak in the low-mass star region that are slightly better fitted 
with a separate BD population and an overlap region. Al- 
though these features alone do not reject the continuous IMF 
model, they might be taken as a further support of the argu- 
mentation towards a two-populations model. 

It should be mentioned that the continuous IMF model (Fig. 
|6l first panel) still fits the observed Trapezium and Pleiades 
MF with high confidence (while failing for the other clusters) 
but leads to VLMS binary properties that are inconsistent with 
the observed properties as mentioned in §[T] A possible exten- 
sion of our modelling would thus be to fit both the observed 
IMF and the observed binary statistics. 

4.2. BD to Star Ratio 

We analyzed the BD-to-star ratio TZ of TA, Trapezium and 
IC 348, defined as 



7^: 



A^(0.02-0.075Mq) 
A?(O.15-1.OM0) 



(14) 



where is the number of bodies in the respective mass 
range. The mass ranges are chosen to match those used in 
iKroupa & Bouvier (2003) in their definition of TZ. Since BDs 
below 0.02 Mq are very difficult to observe and since TA and 
IC 348 do not host many stars above 1 Mq (in contrast to 
Trapezium and the Pleiades), we restricted the mass ranges to 
these limits. 

The Pleiades cluster is difficult to handle due to a lack 
of BD data. Moreover, at an age of about 130 Myr 
jBarrado v Navascues. Stauffer & Javawardhanall2004l) it has 
undergone dynamical evolution, and massive stars have al- 
ready evolved from the main sequence which affects the 
higher mass end of the IMF (Moraux et al. 2004). 

Another relevant quantity is the ratio T^hbl of BD-like to 
starlike objects at the HBL, i.e. the classical BD-star border. 



7^, 



HBL ■ 



Cbd(0.075Mo) 

estar(0.075MQ) ■ 



(15) 



In the classical continuous IMF approach it is one by defini- 
tion, because otherwise the IMF would not be continuous. In 
a two-component IMF its value depends on the shapes of the 
BD-like and the starlike IMF, as well as on the BD-to-star ra- 
tio. The evaluation of equation ( fTsT i yields T^hbl = 0.17 for 
the Trapezium, TZhbl = 0.30 for TA, TZhbl = 0.22 for IC 348, 
and 7?.HBL = 0.3 for the Pleiades. Although 7?.hbl is not an 
input parameter here but is calculated from the two IMFs and 
their relative normalization, 7?,pop, it could easily be used as 
one instead of 7?.pop while calculating the latter from 7?.hbl- 

Figure |9] shows the contour plots of the 5% and 1% sig- 
nificance ranges in the aBD - T^hbl space for Trapezium, TA, 
IC 348, and the Pleiades. The significance values are cal- 
culated from via the incomplete Gamma function as de- 
scribed in § |3] The contours mark the regions outside which 
the hypothesis of a two-component IMF with a single power 




2 component 
random pairing 
continuous 
Lada (2006) 
Kroupa et al. (2003) 
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Fig. 1 1 . — The binary fraction /(nip,i,„) of our best 2-component IMF fit 
for the Trapezium (solid line) with qbd = 0.5 in comparison with the binary 
fraction, /cCHprim). of the best continuous IMF model (dotted line, see also 
top panel in Figure |6). The 2-component model assumes random pairing of 
stars but equal-mass pairing for BD binaries. The continuous model, on the 
other hand, assumes random pairing of companions chosen randomly from 
the canonical IMF (asD = 0.3, cui = 1.3 and = 2.3). This model is merely 
shown for illustrative purposes and has already been discarded (§ [T). The 
thin dashed line refers to the same two-component IMF as the solid line, 
but with random-pairing instead of equal-mass pairing for BDs. The visible 
jumps in /(»ip,i,^) at log(m/MQ) = — 1.15 and -0.7 are due to the truncation 
at the upper and lower mass limit of the BD-like and the starlike population 
with their different binary properties. The global shapes of the continuous 
and the two-component IMF curves are similar but /c('Mprim) is significantly 
higher than /(mp,i„,) in the low-mass star and VLMS regime and outside the 
uncertainties of the data bv iLadaH2006i) . 



law for BDs and a double power law for stars has to be re- 
jected with 95% or 99% confidence. Also shown (by a cross) 
is the standard configuration with a^D = 0.3 and log7?.HBL = 
for the continuous standard IMF. This point is outside both 
levels for all three clusters, and at least for Trapezium and 
TA it is well outside even for arbitrary aso- In other words, 
the corresponding hypothesis of a continuous IMF has to be 
rejected with at least 99% confidence. 

The size of the non-rejection areas can be used for an es- 
timate of the errors of obd and 7?,hbl- However, one has to 
keep in mind that these are only maximum possible deviations 
with all the other parameters kept at the optimum. The non- 
rejection areas in the full three-dimensional parameter space 
are therefore expected to be somewhat narrower 

4.3. Binary Fraction 

Additionally, the binary fraction for each cluster as well as 
the total binary fraction is calculated for the best-fit models. 
The fraction of binaries as a function of the primary mass, 
/(mpi-im), among all systems is shown in Figure [TO] For stars 
the binary fraction is a monotonic function of the primary 
mass in agreement with the data (lLadal l2006), at least for the 
Trapezium, IC 348, and the Pleiades. For BDs, however, it 
is flat due to the equal-mass pairing. In the case of random 
pairing it would approach zero for very low mass BDs (Fig. 
fTTT i. The true mass ratio function grows monotonically with 
th e mass ratio and becomes very steep near ^ = 1, as shown 
bv lReid et al.l (12006) . Thus, for BDs the true binary fraction is 
probably closer to the equal-mass case than the random pair- 
ing case. 

For comparison. Figure [TT] also shows the binary fraction, 
/c('«prim), for the Trapezium that would result from a contin- 
uous IMF. Although the overall shape is very similar to that 
of the two-component model the binary fraction near the BD- 



10 



Thies & Kroupa 



star transition is significantly higher for low-mass stars than 
the observed values while being approximately equal for stars 
above 1 Mq. Thus, a continuous IMF cannot fit the obser- 
vational data as good as a two-population IMF even if /tot is 
reduced. 

The total fractions of BD-BD binaries, /bd-bd, of star-star 
binaries, /star-star, and the fraction of (very low mass) star- 
BD binaries, /star-BD, are of further interest. Pairs of the 
latter type consist of two objects of the BD-like or stailike 
population (see § |5] for the motivation and definition of the 
populations) but where the primary object is a star (mprim > 
O.O75M0) within the BD-VLMS overlap region between 0.07 
and 0. 15 Mq, while the companion is a true, physical BD with 
a mass below 0.075 Mq. 

For each cluster we define 

„ A^BD-BD 

/bd-bd = t;^ (lo) 



A^sys,BD 
A^MD-MD 



sySjMD 



star— star ~ 



-^sys,star 



/sta 



star-BD 



star-BD = 



(17) 



(18) 



(19) 



*sys,star 

where A^bd-bd is the number of all BD-BD binaries (i.e. all 
objects have masses < 0.075 M0), A^md-md that of all M 
dwarf-M dwarf (MD-MD) binaries, A^stai-stai that of all star- 
star binaries, and A^star-BD the number of mixed (star-BD) bi- 
naries (composed of a VLMS and a BD). Furthermore, A^sys.BD 
is the number of all BD systems (including single BDs), 
A^sys.MD that of all M dwarf systems, and A^sys.star the num- 
ber of all systems with a star as primary and with a stellar, 
BD, or no companion. As for the BD-t q-star ratio we also 
applied the (primary) mass ranges from iKroupa & Bouvier 
(12003) on the BD and star sample from each cluster but with 
no gap between the BD and the stellar regime, i.e. O.O2M0 < 
wjprim < 0.075Mo for BDs and O.O75M0 < nip^^ < IMq for 
stars. In addition, the binary fraction of M dwarfs is calcu- 
lated in the same way as the stellar one but in the mass range 
0.075 M0 < mpi-im < 0.5 M0. The uncertainty limits of the bi- 
nary fractions and BD-to-star ratios have been derived from 
those of the IMF slopes by applying the minimum and maxi- 
mum slopes. 

The results are shown in Table |4] The binary fractions for 
BDs vary only slightly between 12% and 15%, while the stel- 
lar binary fraction, /star-star, is about 70% for TA and 30%- 
40% for the others. They are apparently slightly lower than 
the binary fractions that are set for the starlike population be- 
cause of the nonconstant distribution binary as shown in Fig- 
ure [TO] The binary mass function (eq. ||6l) is smaller in the 
mass region below 1 Mq and thus the binarity is below aver- 
age if the focus is set on this region. 

Furthermore, in the overlap region the relatively low binary 
fraction of VLMSs from the BD-like population also con- 
tributes to /star-star, which rcsults in an even lower value of 
/stai-stai - This trend is more emphasized for the M dwarf bi- 
nary fraction, /md-md, which is about 10% lower than /star-star 
for each cluster but still much larger than /bd-bd- 

The star-BD binary fraction /star-BD is of special interest 
since it is the measure for the "dryness" of the BD desert. 
Note that due to the equal-mass pairing used here for BD-like 
binaries, the BD-like population formally does not contribute 
to /stai-BD- All star-BD binaries are from the starlike regime 



which extends down to wo star = 0.07 Mq, i.e. into the BD 
mass regime. For the Trapezium, IC 348, and the Pleiades our 
two-component models yield values between 2% and 2.5%, 
whereas TA shows /star-BD ~ 5%. For comparison, the con- 
tinuous IMF from the top panel of Figure |6] corresponds to 
/star-BD = 15% ±4%. For mass ratio distributions other than 
equal-mass pairing we expect a higher /star-BD, consisting 
mostly of binaries from the BD-like population with primary 
masses slightly above the HBL and companion masses below. 
For random pairing of BD-like binaries this increment is be- 
tween 0.01 and 0.03. Also the lower mass limit of the starlike 
population slightly influences /ot. It should also be noted that 
the size of the error limits for /star-BD are calculated from the 
uncertainties of the IMF model (Tab. |3]l but do not include the 
uncertainties of mo.star- 

5. DISCUSSION: BROWN DWARFS AS A SEPARATE 
POPULATION? 

5.1. An Apparent Discontinuity 

By correcting the observed MFs for unresolved multiple 
systems a discontinuity in the IMF near the BD/VLMS region 
emerges. We have also tried to model continuous single-body 
IMFs, but we find this hypothesis of continuity to be incon- 
sistent with the observed MFs given the observational data on 
the binary properties of stars and BDs. We have shown that 
the empirically determined difference in the binary properties 
between BDsA'LMSs on the one hand side and stars on the 
other, and the empirical finding that stars and BDs rarely pair, 
implies a discontinuity in the IMF near the BD/VLMSs mass. 

Thus, the discontinuity in binary properties, which has 
already been interpreted to mean two separate populations 
dKroupa & Bouvierll2003h . also implies a discontinuity in the 
IMF. This relates the probably different formation mechanism 
more clearly to the observational evidence. 

We have performed a parameter survey allowing the IMF 
parameters ubd, 7?.pop, and mmax.BD to vary finding that the 
canonical stellar IMF (ai = 1.3 and a2 = 2.3) cannot be dis- 
carded even for TA, and that the BD/VLMS discontinuity is 
required for all solutions. The discontinuity uncovered in this 
way, if measured at the classical BD-star border, is of a simi- 
lar magnitude for the stellar clusters studied (0.17 < 7<^hbl < 
0.30), supporting the concept of a universal IMF, which is, by 
itself, rather notable. 

We recommend calling these populations BD-like and star- 
like with respect to their formation history. These populations 
have probably overlapping mass ranges since there is no phys- 
ical reason for the upper mass limit of the BD-like population 
to match the lower mass limit of the starlike one. Furthermore, 
the best-fit models suggest such an overlap between 0.07 Mq 
and about 0.2 Mq. According to this classification through 
the formation history, BD-like objects would include VLMSs, 
while starlike ones would include massive BDs. 

The overlap region implies that BD-like pairs can consist 
of a VLMS-BD pair, and that a starlike binary can consist 
of a stellar primary with a massive BD as a companion. As 
can be seen from Table 21 the star-BD fractions are similar 
(except for the dynamically unevolved TA) and that the star- 
BD binary fraction, /star-BD, is about 2%-3%. This is some- 
what highe r than the value of < 1 % BD companion fraction 
Grether & LineweaveJ (l2006h found for nearby stars but still 
far below the value expected for a single population model (at 
least 7%-15% in ,Ki-oupa et al..2003iand 15% ±4%inourcal- 
culations). Simulations bv lBate etall ( l2003h . lBate & Bonnelll 
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TABLE 4 

Binary fractions /bD-BD for BD-BD, /mD-MD for MD-MD, /stai-BO for STAR-BD and /stai-slar FOR STAR-STAR BINARIES IN THE STUDIED 
CLUSTERS. Note that the STAR-BD binaries result from the PAIRING IN THE OVERLAP REGION BETWEEN ABOUT 0.075 AND O.ISMq. THUS, 

FOR EXAMPLE, IN IC 348 2.1% OF ALL STARS ARE EXPECTED TO HAVE A PHYSICAL BD AS A COMPANION. IN ADDITION, THE BD-TO-STAR RATIO 
BETWEEN 0.02 AND I.OMq TZ, AND AT THE HBL, T^HBL, IS LISTED. NOTE THAT THE NUMERICAL UNCERTAINTIES OF /s,ar-BD ARE PROBABLY MUCH 

SMALLER THAN TRUE ONES SINCE /s,ar-BD LARGELY DEPENDS ON niQ star. 



Cluster 


/bD-BD /mD-MD /star-star /star-BD 


TZ 72,HBL 


Trapezium 
TA 

IC 348 
Pleiades 


0.13±0.01 0.30±0.01 0.34±0.01 0.023±0.002 
0.15 ±0.01 0.64 ±0.06 0.69 ± 0.05 0.046 ±0.005 
0.13 ±0.02 0.29 ±0.02 0.33 ±0.02 0.021 ±0.004 
0.13 ±0.02 0.33 ±0.03 0.37 ±0.03 0.025 ±0.002 


0.18 ±0.03 0.17 ±0.04 
0.27 ±0.09 0.30 ±0.11 
0.20 ±0.10 0.22 ±0.11 
0.28±0.18 0.28±0.18 



(HOOI), and lBatd (l2005h predict a fraction of about 2% star- 
BD binaries (actually one M dwarf with a BD companion out 
of 58 Stars formed in three independent calculations). Note 
that, in our model, /stai-BD is a prediction of the required over- 
lap region and is sensitive to the overlapping range. As our 
model for the Pleiades IMF suggests, the overlapping range 
might be considerably smaller than that we have found for 
the Other clusters. Furthermore, we did not fit the lower mass 
border, mo.star, of the starlike population but simply assumed a 
value of 0.07 Mq which is well in the BD mass regime (and is 
the major source of star-BD binaries in our Pleiades model). 
A slight increase of mo.star by only 0.01 Mq would cause the 
Pleiades star-BD binary fraction to drop to nearly zero. 

Several authors doubt the existence of two separate popu- 
lations. Most recently, Eisloffel & Steinacker (2007) summa- 
rize that the observational community in general prefers the 
model of starlike formation for BDs. They mention the detec- 
tion of isolated proto(sub)stellar "blobs" in the Ophiuchus B 
and D clouds, which may support the theory of starlike for- 
mation for BDs. It remains unclear, though, how many of 
these blobs will actually form BDs in stead of dissolving from 
lack of gravitational binding energy. iGoodwin & WhitworthI 
(l2007h refer to a private communication with A. Nordlund 
stating that the pure turbulence theory predicts about 20,000 
transient cores for every actual pre-stellar core of about 

0.1 Mq. 

We recall that one of the main reasons for the existence of a 
separate population is the semi-major axis di stribution of BD 
binaries (Fig. [T]). But Luhm an et alj (|2007) also mention a 
wide binary BD in Ophiuchus with a separation of approxi- 
mately 300 AU. Indeed, a small number of wide BD binaries 
are known. However, it can be doubted that the occasional 
discovery of a wide BD binary may expand the narrow semi- 
major axis distribution (Fig. [T]) to a starlike one. The striking 
evidence posed by the lack of BD companions to stars is a 
strong indication for two populations. It is usually ignored by 
the community, though. We note that even if there actually 
may be some BDs that formed starlike they are most probably 
a minority. 

There is also the interpretation of the BD desert being a 
"low-<7 desert" rather than an a bsolute mass-dependent drop in 
the companion mass function. iGrether & Lineweaved (l2006h 
find a low-mass companion desert of solar-type primary stars 
between approximately 0.01 and 0.06 Mq. They find this in- 
terval to be dependent on the primary star mass and therefore 
predict M dwarfs to have BD companions, and that M dwarfs 
ought to have a companion desert between a few Jupiter 
masses and the low-mass BD regime. However, this inter- 
pretation does not address the different orbital properties of 



BDs and stars as well as the different q distribution. 

5.2. Implications for the Formation History 

Can the existence of such a discontinuity, i.e. the forma- 
tion of two separate populations, be understood theoretically? 
Although BDs and stars appear to be distinct populations the 
formation of BDs is likely to be connected to star formation. 
Bate et al. (2002, 2003) show that BDs significandy below 
0.07 Mq cannot form in a classical way since the minimum 
mass they need for stellar-type formation would also lead to 
progressive accretion and growth to stellar mass unless they 
are in regions with very low mass infall rates. But such re- 
gions are very rare, because the prevailing densities and tem- 
peratures c annot achieve the required Je ans masses, as also 
stressed by 'Goodwin & WhitworthI (l2007l) . For this reason, 
lAdams & Fa tuzzo (1996) expected BDs to be rare. Only a 
small fraction of the BDs, especially those at the high-mass 
end of the BD regime, may form this way if the surrounding 
gas has been consumed by star formation processes just af- 
ter the proto-BD has reached the Jeans mass. To explain the 
actually higher BD frequency (per star) in recent surveys the 
accretion process has to be terminat ed or impeded somehow 
(iBonnell. Larson & Zinneckeiil2007|) . 

Also the above-mentioned differences in the distribution 
of the semimajor axes between BDs and stars cannot be ex- 
plained by a scaled-down star formation process, because that 
would imply a continuous variation and a much broader semi- 
major axi s distribution for BD s and VLMSs that has not been 
observed (iKroupa et al.l2003h . While binary stars show a very 
broad distribution of their semimajor axes peaking at about 
30 AU the semimajor axes of BDs are distributed around 
about 5 AU with a sharp truncation at 10-15 AU (Fig. [T]i. No 
smooth transition region between both regimes can be recog- 
nized (C lose et al. 2003). In a high angular resolution survey 
iLaw et al.i (i2007) found that the orbital radius distribution of 
binaries with V-K < 6.5 appears to differ significantly from 
that of cooler (and thus lower mass) objects, suggesting a 
sudden change of the number of binaries wider than 10 AU 
at about the M5 spectral type. This is in agreement with our 
finding of a possible BD-like population that extends beyond 
the hydrogen-burning mass limit into the VLMS reg ime. 

In a radial-velocity survey of Chamaeleon I, iJoergensI 
(l2006l) has found evidence for a rather low binary fraction be- 
low 0. 1 AU, while most companions found in that survey orbit 
their primaries within a few AU. For this reason an extreme 
excess of close BD binaries that cannot be resolved by imag- 
ing surveys appears to be unlikely. A larger binary fraction 
than ab out 15% would thus not be plausible. Basri & ReinersI 
( 120061) suggest an upper limit of 26% ± 1 0% for the BD binary 
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fraction based on their own results (1 1%^q q4, for separations 
below 6 AU) and the survey by Close et al. (2003; 15% ±7%, 
for separations greater than 2.6 AU) by simple addition of the 
results. This is nearly consistent with a BD binary fraction 
of 15%, since the survey is neither magnitude- nor volume- 
limited. However, they admit that their value may be over- 
estimated since the objects with separations between 2.6 and 
6 AU are counted twice. We note further that even a total BD 
binary fraction of 25% (/bd-bd = 0.25), although outside the 
error limits of our best-fit models, would only lead to a minor 
change in the fitted IMFs. 

It has been argued jBasri & ReinersI 12006) that the lower 
binary fraction of BDs is just the extension of a natural trend 
from G dwarfs to M dwarfs (Figs. [TO] and \TT[ . Our con- 
tribution has shown that this trend can be understood by the 
simple fact that there are many fewer possibilities to form a 
binary near the lower mass end than for higher component 
masses. The observational data are in better agreement with a 
minimum mass, »Jmin, near the hydrogen-burning mass limit 
and a low overall binary fraction of BD-like objects than with 
an "all-in-one" IMF from the lowest mass BDs to the up- 
per stellar mass limit, as shown in Figure [TT] This observed 
trend thus appears as an additional enforcement of the two- 
populations model of BDs and stars. 

Given that the conditions for a starlike formation of BDs are 
very rare (Bate et al. 2003), four alternative formation scenar- 
ios for BDs apart from starlike formation can be identified, 
namely 

1 . Formation of wide star-BD binaries via fragmentation 
of a proto- or circumstellar disk and subsequent disrup- 
tion by moderately close encounters. 

2. Formation of BDs as unfinished stellar embryos ejected 
from their birth system. 

3. Removal of the accretion envelopes from low-mass pro- 
tostars via photoevaporation. 

4. Removal of the accretion envelopes due to extremely 
close stellar encounters (Price & Podsiadlowski 1995). 

Scenario 4 can be ruled out as the major BD formation 
mechanism because the probability of such close encoun- 
ters, with r equired flyby distances typically below 10 AU 
dKroupa & Bouvier 2003) for efficient disruption of accre- 
tion envelopes (less than a tenth of those proposed by Thies 
et al. [2005] for triggered planet formation), is far too 
low for such a scenario being a significant contribution to 
BD formation. The pho to evaporation model, as studied 
by IWhitworth & Zinneck er (2004), also cannot be the ma- 
jor mechanism of BD formation ( Kroupa & Bouvieiil2003l) . 
It predicts a variation of the IMF with the population num- 
ber and density of the host cluster. In dense starburst clus- 
ters (young globular clusters) with a larger number of 0/B 
stars or even modest clusters such as the ONC with a dozen 
O/B stars compared to TA, a larger fraction of low-mass stars 
would have halted in growth. This would result in a bias to- 
wards M dwarfs, since many of them would be failed K or 
G d warfs. In contrast to this prediction. 'Bric eno et al.l (l2002h 
and lKroupa et al.l (12003b show that the IMFs of TA and ONC 
are very similar in the mass range 0.1-1 Mq, while globular 
cluste rs likewise ha ve a low-mass MF similar to the standard 
form ( lKroupall2001h . 



5.3. Embryo Ejection 

iReipurthI (|2000|) and lReipurth & Clarkd (1200 lb introduced 
the formation of BDs as ejected stellar embryos as the alter- 
native scenario 2. If a forming protostar in a newborn multi- 
ple system is ejected due to dynamical instability its accretion 
process is terminated and the object remains in a protostellar 
state with only a fraction of the mass compared to a fully de- 
veloped star. Since the final mass is physically independent 
of the hydrogen fusion mass limit one would not expect the 
mass range of ejected embryos to be truncated at the HBL 
and thus expect an overlap region between these populations. 
This fully agrees with the requirement of having to introduce 
such an overlap region in order to fit the observed IMFsy, in 

m 

This model gives some hints to understand the low BD bi- 
nary fraction as well as the truncation of the semimajor axis 
distribution of BDs. The decay of a young multiple system of 
three or more stellar embryos typically leads to the ejection 
of single objects but also to the ejection of a small fraction of 
close binaries. In order to survive the ejection, the semi-major 
axis of such a binary must be significantly smaller (by a factor 
of about 3) than the typical orbital separation within the orig- 
inal multiple system. A similar explanation is that the orbital 
velocity of the BD binary components has to be higher than 
the typical ejection velocity in order to keep the interaction 
cross section of the binary with other system members small. 
Indeed, the veloc ity dispersion of BDs in th e embryo-ejection 
model shown in Kroupa & Bouvierl (l2003h is < 2kms"' for 
the majority of the BDs. This is in good agreement with 
the Keplerian orbital velocity of each BD-binary member of 
about 1.5-2kms~' for an equal-mass binary of 0.05-0.08 Mq 
and a « 10 AU. The majority of BD binaries have smaller sep- 
arations, and, consequently, higher orbital velocities and are 
bound tighter. This would set the low binding energy cut near 
^bind = 0.2pc"^Myr^ in Figure|2l 

There have been nur nerical sirnulation s of star formation 
and dy namics, e.g. by iBate et all ( l2003h and lUmbreit et alj 
(12005), in which binaries are produced via ejection that show 
remarkably similar^ properties to the actually observed ones. 
lUmbreit et al.l ( 120051) describe the formation of BDs from de- 
caying triple systems. Their simulations predict a semi-major 
axis distribution between about 0.2 and 8 AU (see Fig. 8 in 
their paper), peaking at 3 AU. This is slightly shifted towards 
closer separations compared with the results by Close et aQ 
(|2003|) but still i n agreement with the observ ational data. In 
contrast to this, Goodwin & WhitworthI ( l2007i) doubt the fre- 
quent formation of close BD binaries via ejection, arguing 
that hydrogen-burning stars which formed via ejection were 
almost always single. 

In further qualitative support of the embryo ejection model, 
iGuieu et"aLl (1200 6) describe a deficit by a factor of 2 of 
BDs near the highest density regions of TA relative to the 
BD abundance in the less dense regions that can possibly 
be explained by dynamical ejection and consequently larger 
velocity dispersion of stellar embryos, i.e. BDs. A star- 
like fragmentation scenario would result in an opposite trend 
since the Jeans mass is smaller for higher densities, thus al- 
lowing gas clumps of lower mass to form (sub)stellar bod- 
ies. Contrary to this, iLuhmanI (120061) did not find any ev- 
idence for a diff erent spatial d istribution of BDs and stars 
in TA. Recently, iKumar & Schmeia ( 2007i) have found that 
substellar objects in both the Trapezium and IC 348 are dis- 
tributed homogeneously within twice the cluster core radii 
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while the stellar populations display a clustered distribution. 
They conclude that these distributions are best explained with 
a higher initial velo city dispersion of BDs, in accordance with 
iKroupa & BouvieJ |2003), supporting the embryo-ejection 
model. 

However, the embryo-ejection model has a challenge to re- 
produce the high fraction of significant disks around BDs that 
have been observ ed in young clusters. Several studies, e.g. 
iNatta et aLl (120041) . reveal a considerable number of BDs with 
an infrared excess that indicates the presence of warm cir- 
cum(sub)stellar material. While these studies do not show 
the actual mass of these disks because a small amount of 
dust in these disks is suf ficient to produce these excesses, 
IScholz. Javawardhana &"W ood (2006) found 25% of BDs in 
TA having disks with radii > 10 AU and significant masses 
(larger than 0.4 Mj). For the remai ning 75% n o disks were 
detectable. In simulations, e.g. by iBate et al] (|2003), such 
large disks survive occasionally, but les s frequently (abou t 
5%) than sugg ested by the observations (IScholz et al.ll2006h . 
However, Scholz et al.l ([2006) point out that their results to 
not rule out the embryo-ejection scenario and admit that this 
mechanism may still be relevant for some BDs. In gen- 
eral, the embryo-ejection model is in agreement with at least 
the exist ence of low-mass circumsubstellar disk s up to about 
10 A U (Bate et al. 200l lUmbreit et all l200l iGiidel etal.1 
l2007h . Thies et al. (2005|) also show that a thin low-mass 
disk can survive near-parabolic prograde coplanar encounters 
above about three Hill radii with respect to the disk-hosting 
BD. This means that a disk with a radius u p to about 5 AU 
(later viscously evolved to about 10 AU, see ' Bate et al.ll2003h 
can survive a flyby of an equal-mass embryo within about 
15 AU while larger disks or widely separated binaries would 
be disrupted. 

Furthermore, iBate et akl (l2003l) suggest from their simula- 
tions that the binary fraction via ejection might by as small as 
about 5% (see also.Whitworth et al. 2007). 

5.4. Disk Fragmentation and Binary Disruption 

The fragmentation of protobinary disks with subsequent 
disruption of a star-BD binary is another promising alternative 
scenario. A disk can fragment during the accretion process if 
it reaches a critical mass above which the disk becomes gravi- 
tationally unstable against small perturbations. Fragmentation 
may be triggered by an external perturbation, i.e. infalling 
gas clumps or a passing neighbor star The latter mechanism 
may also be capable of triggering fragmentation in relatively 
low- mass circumstellar disks resulting in rapid planet forma- 
tion ( Thies et ani2005h . Thus, BD formation via disk frag- 
mentation Ukely plays an important role in the early ages of 
the cluster where the frequency of massive disks is highest 
dBaisch et al. 2001; Whitworth et al. 2007|). 

Following the argumentation of IGoodwin & WhitworthI 
( 12007 ). disk fragmentation may explain the observed distribu- 
tion of BD binary separations at least as well as the embryo- 
ejection model. In addition, it may explain the existence of 
wide star-BD binaries, since a fraction of the initial wide bi- 
naries can survive without being disrupted. Because the likeli- 
hood of disrupting close encounters depends on the mass and 
density of the host cluster, one expects a higher fraction of 
those wide star-BD binaries in smaller and less dense clusters 
and associations like TA. However, further studies and obser- 
vatio ns are needed to test this hypothesis. 

iWhitworth & Stamatellos (2006) show that BDs can actu- 
ally form as widely separated companions to low-mass stars 



with mass w at a sufficiently large disk radius, rdisk, 

rdisk > 150 AU f — j , (20) 

where the disk is cool enough to allow a substellar clump 
to undergo gravitational collapse. For a primary star below 
0.2 Mq this minimum radius therefore becomes less than 
90 AU which is in rema rkable agreement with th e two wider 
VLMS binaries found bv lKonopackv et al.1 (l2007h . Such large 
distances to the primary star allow the formation of BD-BD 
binaries as well as the survival of circumsubstellar disks up 
to about 10-30 AU, depending on the total mass of the pre- 
substellar core and the mass of the primary star. Furthermore, 
this scenario explains the existence of wide star-BD binaries. 
In addition, such a wide star-BD binary can be disrupted by 
moderately close encounters of about 100-200 AU (i.e. a dis- 
tance similar to the star-BD orbital radius), the disruption of 
such systems appears to be likely in contrast to the disrup- 
tion of accretion envelopes as required in the already rejected 
scenario 4. 

5.5. Summary 

Both the embryo-ejection model and disk fragmentation 
with subsequent wide binary disruption explain the above- 
mentioned connection between stars and BDs, since BDs start 
to form like stars before their growth is terminated due to their 
separation from their host system or from lack of surrounding 
material in the outer parts of a circumstellar disk. It is obvious 
that the formation rate of these embryos is proportional to the 
total star formation rate. 

For these reasons these formation mechanisms appear to 
be the most likely ones for BDs and some VLMSs and 
BD/VLMS binaries. It cannot be decided yet which scenario 
is the dominant mechanism. This may depend on the size and 
the density of the star-forming region. We expect, on the other 
hand, the classical starlike formation scenario to be of some 
importance only for the most massive BDs. 

The possibility of two different alternative BD formation 
mechanisms (disk fragmentation and embryo ejection) may 
lead to another discontinuity in the intermediate mass BD 
IMF since both scenarios correspond to different binary frac- 
tions as mentioned above. The currently available data, how- 
ever, are far from being sufficient for a verification of this 
prediction. 

6. CONCLUSIONS 

The different empirical binary properties of BDs and stars 
strongly imply the existence of two separate but mutually re- 
lated populations. We have shown that if the IMF of BDs and 
stars is analyzed under consideration of their binary proper- 
ties then there is a discontinuity in the transition region be- 
tween the substellar and stellar regime that is quite indepen- 
dent of the host cluster The discontinuity in the IMF near 
the HBL is a strong logical implication of the disjunct binary 
properties and suggests splitting up the IMF into two compo- 
nents, the BD-like and the starlike regime. An alternative but 
equivalent description would be to view the stellar IMF as a 
continuo us distribution function ra nging from about 0.07 to 
150 M© dWeidner & Kroup"al l2004'). and a causally connected 
but disjoint distribution of (probably mostly) separated ultra- 
low-mass companions and ejected embryos with masses rang- 
ing from 0.01 Mq to 0. 1-0.2 Mq. While the canonical stellar 
IMF is consistent with the observed stellar MFs at least for 
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the Trapezium, IC 348, and the Pleiades, the sub-stellar IMF 
of at least the Trapezium and TA has a power law index that 
is consistent with the canonical value a^D = 0.3. Within the 
error limits, our analysis does not reject the canonical power 
law indices for BDs and stars for any cluster. 

The discontinuity is often masked in the observational data 
due to a mass overlap of both populations in the BD-VLMS 
region as well as the higher apparent masses of unresolved 
binaries compared to single objects in the observed IMFjys. 
The discontinuity in the number density near the HBL is a 
step of approximately a factor of 3-5 (Table H]). This im- 
plies a general dependency between both populations and is, 
as far as we ca n tell, consistent with the scen ario of disrupted 
wide binaries dGoodwin & Whitworthll2007b as well as with 
the truncated-star scenario (e.g. as an ejected stellar embryo. 



iReipurth & Cl^llool . since the number of unready stars 
is directly correlated with the total amount of star formation 
in the host cluster. Both embryo-ejection and wide binary dis- 
ruption are also consistent with the properties of close binary 
BDs. 

Our results (Tab. [3]i suggest that about one BD is produced 
per 4-6 formed stars. This suggests the necessity of a dis- 
tinct description of BDs and stars as well as the connection 
between these two populations through their formation pro- 
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